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The quark and gluon viscosities are calculated in quark-gluon plasmas to leading orders in the 
coupling constant by including screening. For weakly interaction QCD and QED plasmas dynamical 
screening of transverse interactions and Debye screening of longitudinal interactions controls the 
infrared divergences. For strongly interacting plasmas other screening mechanisms taken from lat- 
tice calculations are employed. By solving the Boltzmann equation for quarks and gluons including 
screening the viscosity is calculated to leading orders in the coupling constant. The leading logarith- 
^^ ] mic order is calculated exactly by a full variational treatment. The next to leading orders are found 

^\. , to be very important for sizable coupling constants as those relevant for the transport properties 

Qs ' relevant for quark-gluon plasmas created in relativistic heavy ion collisions and the early universe. 
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I. INTRODUCTION 

Transport and relaxation properties of quark and gluon 
(QCD) plasmas are important in a number of a different 
contexts. They determine the time that it takes a quark- 
gluon plasma formed in a heavy-ion collision to approach 
equilibrium, and they are of interest in astrophysical sit- 
uations such as the early universe, and possibly neutron 
stars. 

The basic difficulty in calculating transport properties 
of such plasmas, as well as of relativistic clcctron-photon 
(QED) plasmas, is the singular nature of the long-range 
interactions between constituents, which leads to diver- 
gences in scattering cross sections similar to those for 
Rutherford scattering. This makes the problem of fun- 
damental methodological interest, in addition to its pos- 
sible applications. The first approaches to describe the 
transport properties of quark-gluon plasmas employed 
the relaxation time approximation [y||,^ for the colli- 
sion term. This approximation simplifies the collision in- 
tegral enormously and transport coefficients are related 
directly to the relaxation time. The latter is typically 
estimated from a characteristic cross section times the 
density of scatterers. In Refs. [^,H the divergent part of 
the total cross section at small momentum transfers was 
assumed to be screened at momentum transfers less than 
the Debye momentum. However, Debye screening influ- 
ences only the longitudinal (electric) part of the QED and 
QCD interactions, and the transverse (magnetic) part is 
unscreened in the static limit. 

Recently it has been shown that the physics responsible 
for cutting off transverse interactions at small momenta 
is dynamical screening Q. This effect is due to Landau 
damping of the exchanged gluons or photons. Within 
perturbative QCD and QED rigorous analytical calcula- 
tions of transport coefficients to leading order have been 
made for temperatures high [^jj^ as well as low ||] com- 
pared with the chemical potentials of the constituents. 

Transport processes depend on a characteristic relax- 
ation time, Ttr, of the particular transport process con- 
sidered. For example, in high temperature plasma the 
viscosities, r/i = WiTri^i/5, of particle type i are propor- 
tional to the characteristic times for viscous relaxation, 
Trj,i ~ Ttr, which wcrc first calculated in y] to leading 
order in the coupling constant. More generally one finds 
that the typical transport relaxation rates, that deter- 
mines momentum stopping, thermal and viscous relax- 
ation, is in a weakly interacting QCD plasma 
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where the expansion is in terms of the fine structure con- 
stant as = g'^/AiT. The coefficients of proportionality to 
the leading order in ag (in the following called the lead- 
ing logarithmic order) has been calculated analytically 
for a number of transport processes in high temperature 
plasmas [Hja]. Likewise in a QED plasma the typical 



transport relaxation rates for viscous processes, momc 
tum stopping, thermal and electrical conduction have t 
same dependence as (|]) on the QED fine structure cc 
stant a g. 

The dependence of the transport rates on the coupli 
constants is very sensitive to the screening. Besides t 
factor aj. from the matrix element squared of the qua 
and gluon interactions, the very singular QCD inter; 
tions for small momentum transfers lead to a logarith 
^'O-ilrnax / Qmin) , of the maximum and minimum mome 
tum transfers. The typical particle momenta limits t 
maximum momentum transfer, Qmax ~ T, and Det 
and dynamical screening leads to effective screening 1 
small momentum transfers of order qmin ^ <1d ^ 9 
This gives the leading logarithmic order in the coupli 
constant, ln(r/g£)) ~ ln(l/as), to the transport ral 

(0)- 

The calculations in |Q,^ were brief and dealt only wi 

the leading logarithmic order in the coupling consta 

with a given ansatz for the distribution function. He: 

more detailed calculations of the quark and gluon v 

cosities in the high temperature quark-gluon plasmas £ 

presented. The leading logarithmic order is calculat 

exactly by a variational method and the next to leadi 

order - the a1 term in (|l|) - is calculated as well. Becai 

tts is not exponentially small, the next to leading ore 

is important in many realistic physical situations as r 

ativistic heavy ion collisions and the early universe. Fi 

thermore when the Debye screening length is larger th 

the interparticle screening, which occur when aj^O.l 

we shall see below, Debye and dynamical screening brea 

down. Instead lattice gauge calculations have found tt 

quark-gluon plasmas seem to develop a constant scree 

ing mass, nipi ~ I.IT, for temperatures Tx;,2 — 3Tc and 

is important to see what effects this alternative screeni 

mechanism has in strongly interacting plasmas. 

We shall first describe in section II the transport thee 

we use, namely the Boltzmann equation, and the scree 

ing of long range QCD and QED interactions. In secti 

III, we describe the process of shear fiow and the var 

tional calculation necessary in order to find the viscosi 

In section IV we then evaluate the collision term to lee 

ing logarithmic order with a simplifying ansatz for t 

trial function and refer to Appendix A for a full and e 

act variational calculation. In section V we calculate t 

viscosity to higher orders in the coupling constant a 

discuss strongly interacting plasmas. Finally, in secti 

VI we give a summary and discuss generalizations of t 

methods developed here to other transport coefficienti 



II. TRANSPORT THEORY 

Transport processes are most easily described by t 
Boltzmann equation 
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X [nin2(l ± n3)(l ± 724) - (1 ± ni)(l ± 712)713714] 
X S{ei + £2 - £3 " e4)<5pi+p2;p3+P4 7 (2) 

where e, is the energy and pi the momentum of the quasi- 
particles, F some force acting on the quasi-particles, and 
the right hand side of ^ is the colhsion term. 7ii(pi) are 
the Fermi and Bose quasi-particle distribution functions 
for quarks and gluons and the signs ± include stimu- 
lated emission and PauH blocking. The spin and color 
statistical factor 7^2 is 16 for gluons and 12Nf for quarks 
and antiquarks with Nf flavors. |A/p is the squared ma- 
trix element for the scattering process 12 -^ 34, summed 
over final states and averaged over initial states. It is 
related to the Lorentz-invariant matrix element \Ai\'^ by 
|Mp = |A^|V(16eie2e3e4). For gluon-gluon scattering 
§ (see Fig. (0)) 
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where s, t, and u are the usual Mandelstam variables. 
In Eq. (^ the double counting of final states has been 
corrected for by inserting a factor 1/2. For quark-gluon 
scattering 



\M,,f=g\u' + s')[-- — 
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and for scattering of two different quark flavors 
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The matrix element for scattering of the same quark fla- 
vors or quark-antiquark scattering is different at large 
momentum transfer but the same as (pi) at small mo- 
mentum transfers. 

The i~2 aj^fj ^-2 singularities in Eq. (§j|) lead to 
diverging transport cross sections and therefore vanishing 
transport coefficients. Including screening, it was shown 
in [Q,^!^] that finite transport coefficients are obtained. 
In fact, the leading contribution to transport coefficients 
comes from these singularities. In the t — uj^ — q^ channel 
the singularity occurs for small momentum q and energy 
to transfers (see Fig. (|l])). 

For small momentum transfer, q <^ ei, €2 ~ T, energy 
conservation implies that w = ei — €3 ~ vi • q = — V2 • q 
where v^ = Pi. Therefore the velocity projections trans- 
verse to q have lengths |vi_r| = |v2,t| — \/l — fJ,'^, where 
IJ, = oj/q. Consequently Vi^t-V2,t = (1 — ^t^^cos^, where 
(j) is the angle between Vi_t and V2,t- For g <C T we thus 
have 



s ~ -u ~ 2pip2(l ~ cos 6*12) 

~ 2piP2{l — ^)(1 — COS( 
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and the interactions splits into longitudinal and trans- 
verse ones, M 
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The interactions are modified by inclusion of the glue 
or photon, self energies, Hl and IIt [|| (see also Fig. ([ 
In the random-phase approximation the polarizations £ 
given in the long wavelength limit {q <^ T) by 
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where fi — uj/q and ^d — 1/Ad is the Debye wavenumb 
In a weakly-interacting high temperature QCD plasn 

ii 

ql=4TTil + Nf/6)a,T\ (] 

where a^ = g'^ /4:Tr is the fine structure constant for stro 
interactions, the factor {l + Nf/6) is the sum of contrit 
tions from gluon screening, the "1," and from light m; 
quarks, of number of flavors, Nf. In a high temperati 
QED plasma, q^ = inaT^/S, where a is the QED fi 
structure constant. 

One should keep in mind that the self energies of ( 
and (P) are only valid in the long wavelength limit, i. 
for q <C T. When q ^ T other contributions of ore 
OsQT enter (see, e.g., |^) which may be gauge depende 
P4| . However, as long as a^ is small all contributic 
from the self energies can be ignored in the gluon pre 
agator when q ~ T because the matrix element squar 
already carry the order aj.. 

In the above derivations we have consistently assum 
that the screening was provided in RPA by the glu 
self energies which give the Debye and dynamical scree 
ing of longitudinal and transverse interactions respi 
tively. Both effects provide a natural effective cutoff 
momentum transfer less than qmin ~ Qd- These perti 
bative ideas must, however, break down when the scree 
ing length becomes as short as the interparticle sp; 
ing, i.e. q^ ~ T, or in terms of the coupling consta 
a,<{4:n{l+Nf/6))-^ ~ 0.1 according to Eq. {^. InL 
tice gauge calculations of quark-gluon plasma above tw: 
the temperature of the phase transition, Tc ~ 180Me 
one finds strong nonperturbative effects in the plasma 
that the typical screening mass is 7TT,p; ~ I.IT |12]. O 
may argue (l3| that perturbation theory still applies J 
large momentum transfers so that the matrix elemei 
are given by the simple Feynman tree diagrams, but tt 
perturbation theory does not apply for small momenti 
transfers of order q^qo and that one should rather ins( 
the effective cutoff found by lattice gauge calculations 



Ul ~IIt ~ rripi , 
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The phenomenological screening mass of (|11|) provides 
with a method to extend our calculations of transport e 
efficients to larger values for as and it can be combined 
the Debye and dynamical screening in weakly intcracti 
quark-gluon plasmas. 



III. THE VISCOSITY 

With screening included in the interaction we can now 
proceed to calculate transport properties as the viscosity. 
In the presence of a small shear flow, u{y), in the x- 
direction we obtain from the Boltzmann equation 

dni dux sr^ 2 
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X [nin2{l ± rt3)(l ± "-4) - (1 ± '^l)(l ± "2)"-3«4] 
X 6{ei + £2 - £3 - e4)'5pi+p2;p3+p4 . (12) 

For small u we can furthermore linearize the quasi- 
particle distribution function 



n, = nf ^ + 



dn ^ dux 



der, 
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dy 
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where the local equilibrium distribution function is 

nf ^ = (exp[(e. - u • p,)/T] T 1)"\ (14) 

and $i is an unknown function that represents the devi- 
ations from local equilibrium. By symmetry $ has to be 
on the form 



<^=P,Pyf{p/T), 



(15) 



where now the function / must be determined from the 
Boltzmann equation. Inserting (O in the Boltzmann 
equation we find 
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-(l±ni)(l ±n2)n^ni\ 

X 5{ei + £2 — £3 ^ £4)'5pi+p2;p3+P4 

X ($1 +$2 -$3-$4). (16) 

It is very convenient to define a scalar product of two real 
functions by: 



(^ilV'2) 



-;^2X!^i(p)^2(p)^ 
p '^ 
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Thus Eq. (|l^ may be written on the form \X) — I\<^) 
where \X) — PxVy and / is the integral operator acting 
on $. The viscosity is given in terms of <I> [1^] and can 
now be written as 



V = -^2y"PxVy^^p = (X|$) 
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Equivalently, the viscosity is given from dlq as 

(X|$)2 
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Since (.|.) defines an inner product, the quantity 
(Xl*)2/(vl/|/|\l/) is minimal for * = $ with the minimal 



value r]. Equation (19 is therefore convenient for var 
tional treatment, which will be carried out in Appenc 
A. 

To find the viscosity we must solve the integral eqi 
tion (no for which we have to evaluate 

($1 / 1$) = 27ri/2 ^ |Af|2nin2(l±«3)(l±^4) 
Pi,P2,q 

($1 + $2 - <&3 - *4)%. , . ,, 

-x- ;; ^5(£i + £2 - £3 - £4) • (S 



Momentum conservation requires that p3 = pi + q a 
P4 = P2 — q where q is the momentum transfer. Int] 
ducing an auxilliary integral over energy transfers, w, t 
delta-function in energy can be written 

<5(£i + £2 - £3 - £4) = / duj (5(cos 6*1 - ^ - ) 

J Piq 2pig 



X -^5{cos 02-^1+ ) 

P2q 2^2(7 
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where di is the polar angle between q and pi and 62 
the corresponding one between q and p2 (see Fig. ! 
Consequently, we find 



($1/1$) = 



1 



Att^T 



dq 



duj 



/O J-q 

X / dpip\ni{pi){l±ni{pi+uj)) 

'iq-uj)/2 

X / dp2pln2{p2){i±n2{p2-i^j) 

J(q+u>)/2 

r2-K 
/o 27r 



X / ,r-|Af|^($i +$2 -*3 - ^4)^- 
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This integral equation for <I> has been solved in a fi 
cases under simplifying circumstances. F.ex., in Fer 
liquids the sharp Fermi surface restricts all particle n 
menta near the Fermi surface and with a simplified fo] 
for the scattering matrix element techniques have be 
developed to calculate a number of transport coefficiei 
exactly [Q . For the QCD and QED plasmas the very s: 
gular interaction can, once screened, be exploited since 
allows an expansion at small momentum transfers. Th 
an analytical calculation of the transport coefficients c 
be carried out at least to leading logarithmic order in t 
couphng constant |Q,||,^. 



IV. VISCOSITY TO LEADING LOGARITHMIC 
ORDER 

In Ref. Q , through solution of the Boltzmann kine 
equation, the first viscosity of a quark-gluon plasma ■« 
derived to leading logarithmic order in the QCD coupli 
strength. We will in the following give a more thorou 
and exact derivation of the quark and gluon viscosi 



The total viscosity, to leading order, is an additive sum 
of the gluon and quark viscosities, rj = r]g + rjq. 

The leading logarithmic order comes from small mo- 
mentum transfers because of the very singular matrix 
element {uh dominates. For small q the kinematics sim- 
plify enormously and, as we will now show, the integrals 
separate allowing almost analytical calculations. First, 
we can set the lower limits on the pi and p2 integrals 
to zero, however, then replacing the upper limit on q by 
the natural cutoff from the distribution functions which 
is qmax ^ T . Thus wc find from (|22 
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(^1^1$) - ^^8^ y dpiplnx{l±m) 
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X / dp2P2n2{l±n2) 

X |M|2($i+$2-$3-$4)'. 



(23) 



to leading logarithmic order 

The solution to the integral equation or equivalently 
the variational calculation of dl3) is quite technical and 
is for that reason given in Appendix A. A much simpler 
calculation is to make the standard assumption in viscous 
processes, i.e., to take the trial function as 



f{p/T) = (p/Tf 



(24) 



As will be shown in Appendix A this turns out to be a 
very good approximation. It is accurate to more than 

99% for reasons also explained in the appendix^ f can 

be defined up to any constant which cancels in (|l9|) and 
therefore never enters in the viscosity. 

The quantity (<i>i -|- $2 — $3 — $4)^ can be averaged 
over X- and j/-directions while keeping fi and (/) fixed. 
This corresponds to keeping the relative positions of the 
three vectors q, p, and p' fixed relative to each other and 
rotating this system over the three Euler angles (see also 
Appendix A). Consequently, we obtain 



((*i+$2 -$3-$4)^) 



Y^[3(P2-pi)^ 

^(q-(p2-Pi))'] 



r[(3 + /?)(p2 



-pD 



15T4' 

-2pip2(V + 3(1 -Ai^) cos (/))]. (25) 



The integrals over pi and p2 in ([2^) are elementary. Next 
we perform the integrations or averages over fj, and d> 
required in ([2^). We note in passing that the term in ( ps]) 
proportional to piP2 vanishes and that fi^ effectively can 
be replaced by 1/3 (see Appendix A). Let us first consider 
the case of gluon-gluon scattering inserting | AIgg \ ^ from 
(13). We thus find 



mm) 



2V 2 , ri. 
45 ^ 



q dq / dfj, 
'0 

1/2 



Iq^+ULipW \q'+TlTi^i)/{l~^^^W 



This integral is discussed in detail in appendix B. I 
the longitudinal interactions Hl — q^) due to Dot 
screening and the leading term is a logarithm of t 
ratio of maximum to minimum momentum transf 
^T^iqmax / qmin) '^ ^^{T/qjj). Likewise for the transvej 
interactions n^ — J(7r/4)/ig|j due to Landau darr 
ing and the dependence on /x = uj/q provides sufficie 
screening to render the integral finite and the leadi 
term is the same logarithm as for the longitudinal inti 
actions. Whereas the details of the screening are uni: 
portant for the leading logarithmic order, they are i: 
portant for the higher orders and they are calculated 
detail in appendix B. The final result is thus to leadi 
logarithmic order 



($|J|$) = 



7^3 



2V_ 
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'\n{T/qD)T^ 
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Since ($|X) = (64^(5)/7r2)T3, wc find from (|T|) and (^ 
2^15^(5)2 T3 



Vgg = 
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TT^ a^ln{T/qD) 
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to leading logarithmic order in as — g'^ /An. 

To obtain the full gluon viscosity we must add scatti 
ing on quarks and antiquarks which is calculated ana 
gously and only has a few factors different. Firstly, frc 
(H) and (Q) we see that the matrix element squared 
a factor 4/9 smaller. Secondly, the statistical factor 
1^2 = 12A^/ instead of 16. Thirdly, in integrating o^ 
the factor (pf + p\) in Eq. (psl) we note that the dist 
bution function, 71,2, in Eq. (|23|) is now a fermion 01 
Consequently, the pi and p2 integrations give a fact 
(1/2 4- 7/8)/2 less for gluon-quark collisions as compar 
to gluon-gluon collisions and wc find 



Tlg^iVgg+rigq) ' 



Vgg 



l + lliV//48' 
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In [Q the shghtly different result rjg = rigg/{l + Nf/ 
was obtained. 

The quark viscosity can be obtained analogously 
the gluon one. The quark viscosity due to collisions 
quarks only, rjqq, deviates from rjgg by a factor (4/9)^ 
the matrix elements and differences in having Fermi a 
Bose integrals. By comparing to ( [l9| , |l8[|2(j ) we find 



^99 ~ %9 



(15/16)^ 



2q6 



5^3 



(4/9)2(7/8)(l/2) 



Vgg 



28? 
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Note that the statistical factors v cancel in r]gg and 77^5. 
Including quark scatterings on gluons lead to similar fac- 
tors in in (<I>|/|$), namely a factor (9/4) from the matrix 
element, a factor I^JVIN f from statistics, and a factor 
(8/7 + 2)/2 from Bose instead of Fermi integrals. Thus 



'/<? = 



"^m 



1 + 33/7iV/ 



2.2 



1 + \\Nfl\% 
1 + 77V/ /33 



^/% , (31) 



which for N^ = 2 results in rjq — 4:Arig, a quark viscos- 
ity that is larger than the gluon one partly because the 
gluons generally interact stronger than the quarks and 
partly because of differences between Bose and Fermi dis- 
tribution functions. 



V. VISCOSITY TO HIGHER ORDERS IN as 

The leading logarithmic order dominates at extremely 
high temperatures, where the running coupling constant 
is small, but it is insufficient at lower temperatures. The 
next to leading order correction to the viscous rate in 
the coupling constant is of order a'^. It may be signifi- 
cant because the leading logarithm is a slowly increasing 
function. In the derivation of the leading logarithmic 
order, Eq. (P^); ^^ have been very cavalier with any 
factors entering in the logarithm, which are of order a^ . 
It was only argued that the leading logarithmic order 
^ri{qmax/qrmn) ^ \tl[T / qo) bccausc qynax and qrarn wcrc 
of order ^ T and ^^ qjj respectively. Finally, if thermal 
quark-gluon plasmas are created in relativistic heavy ion 
collisions at CERN and RHIC energies, the temperatures 
achieved will probably be below a GeV. We can thus 
estimate the interaction strength from the running cou- 
pling constant as ~ 67r/(33 — INf) ln(r/A) which, with 
A ~ ISOMeV and r< IGeV, gives a^ >0.4. For such large 
coupling constants Debye and dynamical screening is re- 
placed by an effective screening mass, rupi, as discussed 
above which will affect the viscosity considerably. 

To calculate the visco sity to order al exactly, the 5- 
dimensional integral of (|22| ) must be evaluated numeri- 
cally and at the same time a variational calculation of $ 
must be performed. This is a very difhcult task and we 
shall instead use the information obtained in the previ- 
ous section, that the trial function / oc p^ is expected to 
be an extremely good approximation. With that ansatz 
for the trial function, it is then straight forward to calcu- 
late the integral of (p^) numerically and find the viscosity 
to order a^ for the given screening mechanism. The 5- 
dimensional numerical evaluation of the collision integral 
of ( pi| ) is a complicated function of the couphng constant. 
It is convenient to write it in terms of the function Q 
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where the index gg refers to gluon-gluon scattering but 
the analogous definitions applies to gluon-quark and 



quark-quark scattering. The function Q and the effi 
five maximum and minimum momentum transfer, q„ 
and qmin, are given in Appendix B. In weakly intera^ 
ing plasmas, where the screening is provided by Det 
and dynamical screening, the function Q is basically ji 
a logarithm of the ratio of the maximum and minimi 
momentum transfer, i.e.. 
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, a,<0.1 
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By numerical integration we find that the distributi 
functions leads to an effective cutoff of qmax ~ 3r. T] 
is because the distribution functions are weighted wi 
several powers of particle momenta and thus contribi 
the most for p ~ 3T. The effective cutoff is slightly lar^ 
for quark-gluon and quark-quark scattering because t 
Fermi distribution functions emphasize large momer 
than the Bose ones. Debye and dynamic al sc reening lea 
to qmin — 1.2651) as described in Eq. (B6) and so frc 
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0.44 
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The numerical factor inside the logarithm, which gi"\ 
the order a^, is discussed in more detail in Appendix 
In the other limit, qo^T or equivalently a^^O.l, pi 
turbative ideas breaks down and we assume an effecti 
screening mass taken from lattice calculations, qmin 
I.IT, as described by Eq. (|ll|). Thus we find (see (pl( 



Q{ctt/qr,^^n) = 0.626 , 



LO.l. 
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and similarly for quark-gluon and gluon-gluon scatti 

ing Qiqmal/qmzn) = 0.819 and Qiqirial/qmin) = 1-0 
respectively. 

Adding gluon-gluon and gluon-quark scatterings 
obtain the gluon viscosities 
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which extends Eq. 
teracting plasmas 

77g~ 0.342^ [In 



to higher orders. In weakly 
reduces to 
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as{l + Nfl&) 
0.72 



a,(l + iV//6) 
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to leading orders in g^. In s trongly interacting plasn: 
we obtain by inserting (p310|) in (p2|) 



T3 
r]g ~ 0.55 -^ 



1 + 1.31- 



llA^f 
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, a,>0.1 
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In Fig. (^ we show the gluon viscosity with the various 
assumptions for screening. With dash-dotted curve the 
result of Eq. ( pq ) assuming a constant screening mass, 
rUpi = I.IT, is shown. With dashed curve the numerical 
result assuming Debye and dynamical screening of Eqs. 

tand ^ is shown. For a^^O.OS it is given by Eq. 
) to a good ap proximation whereas for as<L0.05 the 
result of Eq. (B4) is better. The final viscosity shown 
by full curve is obtained by combining the two limits, 
i.e., applying Debye and dynamical screening in weakly 
interacting plasmas when qd^T or equivalently a^^O.l 
but an effective screening mass rUpi — I. IT as given by 
Eq. ( pi] ) when as^O.l. This corresponds to chosing the 
smallest value of the viscosities as seen in Fig. (||), i.e., 
the two limits of Eqs. @ and (||). 

Similarly, adding quark-quark and quark-gluon scat- 
terings we find the quark viscosity 
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-iVf- 



[Q 



'Jag) ' 

Hrnax 

Id 



33 



'im) \ 

Hmax I -| — 1 

(Id I 



(39) 



which in weakly interacting plasmas gives 

0.72 



??, ~ 0.7527V/^ [In 



al' \as{l + Nf/Q) 



7Nf 



In 



l.f5 



33 "'\asil + Nf/6) 
and in the strongly interacting plasmas 



,a,<0.1, (40) 



T]g~0.92Nf — 



1 + 1.25!^ 
33 



,a,>0.1 



(41) 



The quark viscosity increases with the number of quark 
flavors, Nf, whereas the gluon viscosity decreases as can 
be seen in Fig. (Q) and (^, where the viscosities are 
shown for two and three flavors respectively. The total 
viscosity of a quark-gluon plasma, r/ = r/g + rjq, is domi- 
nated by the quark viscosity. 

From the definition of the viscosity in terms of the col- 
lision integral (Iq) and ( pO| ) , which only contains positive 
quantities, it follows trivially that the viscosity is positive 
as is a physical necessity. The resulting viscosities of Eqs. 
( |3q ) and ( |39| ) are positive quantities whereas the ag^O-l 
expansions of Eqs. ( pq ) and (Efl) are not when extended 
to the region a^^O.S. This explains the results found in 
PH l , where it was claimed that estimates of the next to 
leading order a^ could lead to a negative viscosity. 

Contributions from vertex corrections should also be 
considered. In fact for the calculation of the quasi- 
particle damping rate, 7p, Braaten and Pisarski ||I^ 
found that vertex corrections contributed to leading or- 
der 7p?,o — 6.6as for zero gluon momenta, p. Vertex 
corrections do also contribute to order a^ for large quasi- 
particle momenta, p » gT, but they can here be ignored 



since the leading order is jp = Sag ln(l/as) as explain 
in |16|. For the viscosity, however, vertex corrections c 
be ignored since the extra vertices adds a factor a^. Ev 
though integration over soft momenta may cancel a it 
tor as the result is still of higher order in the coupli 
constant. 

Writing each of the viscosities rji {i = q, g) in terms 
the viscous relaxation time, r^i, as 



(327rV45)T4 and 



(4 



where Wg = {32tt^ / A5)T'^ and Wg = {Nf7TT^/15)T'^ i 
the gluon and quark enthalpies respectively, we obtc 
the viscous relaxation rate for gluons 



1 



^9 T 
3353^(5)2 ^ 



Q 



' (99) ' 

Hniax 



llNf 



' (99) " 

Hmax 



alT, 



(4 



and quarks and antiquarks 
1 Il7r9237 T 



3755^(5)2 al 



Q 



' {9q) ' 

HTTiax 



Ilk, 
33 



'ill) ' 

^max 



alT. 



(^ 



The viscous relaxation times, r^^g, Trj^q and t,, = 1I{t^^^ 
T~q ) are thus very similar to the corresponding viscosit 
when divided by a factor of T*. The curves on Fii 
(Q) and (^ therefore applies to the viscous relaxati 
times (times temperature) as well when divided a fact 
of ~ 1.4 and ~ 0.92A^/ for gluons and quarks respectivt 
according to Eq. ( p2| ) . 

In weakly interacting plasma the viscous rates can 
approximated by 



~4.11a^[ln 



0.44 



as{l + Nf/6) 
f 0.72 
48 ""Vas(l + A^//6) 



lliVf , 



a^3.1 (4 



and 



^ 1.27 a^m^ '•'' 



'■ri,q 



H ^In 



a3(l + iV//6) 
1.15 



33 Va4l + ^//6) 
to leading orders in ots- 



l,a,<0.1, (4 



VI. SUMMARY 

By solving the Boltzmann equation for quarks and g. 
ons the viscosities in quark-gluon plasmas were calc 
lated to leading orders in the coupling constant. Inclusi 



of dynamical screening of transverse interactions, which 
controls the infrared divergences in QED and QCD, is 
essential for obtaining finite transport coefficients in the 
weakly interacting plasmas. The solution of the trans- 
port process was extended to strongly interacting plas- 
mas by assuming an effective screening mass of order 
rupi — I.IT, as found in lattice calculations, when the 
Debye screening length became larger than the interpar- 
ticle distance or when as^O.l. The Boltzmann equation 
was solved exactly to leading logarithmic order numeri- 
cally but the result only differed by less than a percent 
from an analytical result obtained by a simple ansatz for 
the deviation from local equilibrium, $ ex PxPy The 
next to leading orders was also calculated and found to 
be very important for the transport properties relevant 
for quark-gluon plasmas created in relativistic heavy ion 
collisions and the early universe. For as^O.l we find 
Vi = Ci^iT^ /al\n{Ci,2/as) whereas for ag^tO-l we find 



Ci,3T /aj with coefficients C^.j given above. 



The viscosity in degenerate plasmas of quarks, i.e., for 
T <C /Xq was calculated in Q. Several differences were 
found. In the high temperature quark-gluon plasma the 
chemical potential can be ignored and the transport pro- 
cesses depend on two momentum scales only, namely T 
and qo ^ gT. In degenerate quark matter three mo- 
mentum scales enter, namely fig, T, and qu ~ g^J'q, and 
the transport process depends considerably on which of 
qo and T is the larger. In fact for T <^ q^ transverse 
interactions turn out to be dominant in contrast to the 
high temperature quark-gluon plasma where transverse 
and longitudinal interactions contribute by similar mag- 
nitude. Furthermore, the existence of a relative sharp 
Fermi surface allows an almost analytical calculation of 
both the leading (logarithmic) order as well as the next 
order a^. 

The techniques for calculating the viscosities to leading 
orders in the coupling constants can be applied to other 
transport coefficients as well. The leading logarithmic 
orders to momentum stopping, electrical conductivities 
and thermal dissipation in QCD and QED plasmas have 
been estimated with simple ansatze for the distribution 
functions in |q|. Based on the experience with the vis- 
cosity studied here, we do not expect the leading loga- 
rithmic order for these transport coefficients to decrease 
by much when a full variational calculation is performed. 
The next to leading logarithmic order for these trans- 
port coefficients can also be estimated with the experi- 
ence obtained above for the viscosity. Good estimates 
are obtained if one in (B5) replaces qmax by the average 
particle momenta entering the collision integral for the 
relevant transport process and qmin ^y ^ qo- 

A few transport coefficients are, however, different. 
The second viscosity C, is zero for a gas of massless rela- 
tivistic particles |l|] and one cannot define a thermal con- 
ductivity in a plasma of zero baryon number. One can, 
however, consider thermal dissipation processes [g| where 
the leading orders also can be calculated with the above 
methods. The effective soft cutoff will, however, be differ- 



ent for thermal dissipation processes as described in 
because the transport of energy introduce dependent 
on Lo which also is present in the transverse screenii 
nT(w/q). 

All the transport processes discussed above depe 
only on momentum scales from the typical partii 
momentum, qmax ^ T down to the Debye screeni 
wavenumber qmm ^ qo ^ gT which also is the mome 
tum scale for dynamical screening. There is, howev 
a shorter momentum scale of order the magnetic ma 
TTT-mag ^ g^T, at which perturbativc ideas of the quai 
gluon plasma fails ||I^. As shown in IQ the quark a 
gluon quasiparticle decay rates depend on this infrar 
cutoff, rrimag- Furthermore, recent studies iQ find th 
the color diffusion and conductivity also depend on tl 
cutoff and therefore the rate of color relaxation is a fact 
l/us larger than Eq. (|l|). 
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APPENDIX A: EXACT VARIATIONAL 

CALCULATION TO LEADING LOGARITHMK 

ORDER 

In this appendix we solve the Boltzmann equation a 
find the deviation from local equilibrium, $, by a var 
tional treatment of Eq. (18). 

For a general function $ = PxPyf{p/T) we have 



$1 + $2 - $3 - $4 = Pl,xPl,yf{Pl) + P2,xP2,yf{P2) 
-P3,xP3,yfiP3) - P4.xP4,yf{P4) 
= -(ctxPy + ClyPx)fip) - HPxPyhip) 

+ {cLxP'y + q.yPx)f{p') + ^^PxP'yfl{p') , (A 

where we have changed notation to p = pi + q/2 
P3 ~ q/2 and p' = P2 — q/2 = P4 + q/2. We have us 
that the energy conserving (5-functions of (|2l|) impl 
pq — p'p — ji. Furthermore, have defined the functioi 



hip)^p'dif/p')/dp=pr-2f, 



(A 



that vanishes when / oc p^ which was the case for t 
ansatz used in section V. 

For small momentum transfer the matrix element ( 
depends only on energy and momentum transfer and t 
azimuthal angle c/). The ^-functions of (|2l| ) taking care 
energy conservation fixes the polar angles 9i and O2 wi 
respect to q. Thus all the angular integrals for fixed 
and (f> reduces to rotating the three vectors q, p and 



over all Euler angles keeping them fixed relatively to each 
other. Only (<l>i + $2 — $3 — $4)^ depends on the Euler 
angles and the integration or averaging over the three 
Eulerian angles, while keeping the relative positions of 
the vectors q, p and p' fixed, i.e. keeping fi and (f> fixed, 
gives 

3 + ^2 

{cixPy + ciyp^y = (q^p; + ciyp'J^ = _,^ , (A3) 



15 



iPxPy)^ = {p'xP'y)^ = Yg , 

2 

{cicPy + q.yPx){PxPyfJ-) ^ J^fJ-'^ , 



(A4) 
(A5) 

(A6) 



2 

{q^Py + ciyp^){ci^p'y + ciyp'J = jg (3pp' + z"^), (A7) 



(cixPy + ciyPx){%p'J = (q^p; + q;p^)(qyP:.) 



15 



(3pp' - 1) 



(qyP.)(p;p>2) = ^(3(pp')'-l): 



(A8) 



(A9) 



where p — p/cp and p' — p'/^p- Since we assume that 
the plasma temperature is much larger than any of the 
particle masses, the particles are relativistic and p, p' 
and q are unit vectors. The vector product of p and p is 
most useful in terms of ^ and (j) (see Fig. (||) 



PP 



I? + {I- ^^)cos0. 



(AlO) 



Next we integrate over /x and (p. The /i integration 
averages /x^ to 1/3 whereas the (j) integration is weighted 
by a factor (1 — cos0 )^ from the matrix elements. Thus 
we find that (A7-A9) vanishes whereby all combinations 
mixing pi and p2 very conveniently disappear. After av- 
eraging over both Euler angles and /i and <j) we obtain 



(($1 + $2-^3 



$4n-\-(io/2 + /2+4//i) 



p2 15 

g22 
p2 5 



,Af + -yn- (All) 



Let us first consider the pure gluon plasma for which 
12^) gives 



($|/|$) = ^/ln(T/to)r3 
15 



f + -^p'f" 







dn 
dep 



dp, (A12) 



where n = (exp(p/T) — 1) ^ is the gluon distribution 
function. Since 



mx) = 



157r2 



fp' 



dn 



dp, 



(Al 



we find from (nSh 

15 T^ 

77: 



fn X dx 



Sng^ HT/qo) 



1 



(Al 



ir + ^x^nndx 





where x = p/T. As mentioned above, the function f{x) 
determined by minimizing (A14). A functional variati 
with respect to / results in a second order inhomogenec 
differential equation for / 



In" fi 



(Al 



where n"/n' = — (l+2n). C is an arbitrary constant th; 
by rescaling /, can be chosen as C = 2 for convenicnc 
For a: 3> 1 we c an approximate n ~ and so we find t 
solution to ( Alf: ), that docs not increase exponentia 
for x — > 00, to be 

„2 



/(^) 



a; > 1. 



(Al 



For a: <C 1 we can approximate n ~ 1/x and the soluti 
to ( A15| ) that is finite at the origin is 

f(x)^x^(C--liix), 0<a;<l, (Al 

5 

where C ~ 0.7 is a constant tha t can only be determin 
by fin ding the full solution to ( A15 ) and matching it 
(A17). This is done by a numerical Runge-Kutta in 
gration and the result is shown in Fig. (|6|). The viscos: 
is now found by inserting / in ( A14). The exact value 1 
77 thus obtained is only 0.523% less than the approxinij 
value, rjgg, of Eq. (pSl). Since the exact value is a va 
ational minimum, it has to be smaller than that of (E 
is only slightly less because / ~ a^ for large as just 
the ansatz of (E4h and / is mainly sampled over values 
X == p/T ^ 1 because the integrals over pi and p2 in (^ 



have powers 



pUo 



p times np{l 



In Ref. 



a variational calculation with trial functions /(p) oc 
lead to a minimal viscosity for i' = 2 .104. This result 
close to the quadratic power of (AM) but tends sligln 
towards the asymptotic form of (A17) (see also Fig. (I 



It has almost the same slope and curvature as the ex£ 
solution around p = 5T (note that the absolute value 
unimportant since it cancels in the viscosity) . The cor: 
sponding viscosity was 0.364% smaller than that of (& 
i.e., in between the exact result and the ansatz f oc x' 
The above analysis was restricted to a pure gl 
plasma. As mentioned above the distribution functic 
are weighted with several powers of momentum and 
do not find much difference between fermions and boso: 
Therefore the deviation from local equilibrium for quai 
will not be much different from gluons and we can be cc 
fident that the ansatz, <I> oc PxPy, of Eq. (|2J) will be 
good approximation for quarks as well accurate witl 
less than a percent. 



APPENDIX B: SOFT AND HARD 
CONTRIBUTIONS 

The essential contribution to ($|/|<i>) is the integral 



Q 



Qmaa 



1 l-q^r 

-1 Jo 



1/2 



nL(M)P 



\q^ + nT{p)/{l-^l^W' 



(Bl) 



For dimensional reasons the function Q can only depend 
on the ratio of Qmax to the momentum scale, qmin, which 
is provided by the screening. For Debye and dynami- 
cal screening qmm ~ Qd whereas lattice calculations of 
strongly interacting plasmas give qmin ^ TTipi = I-IT. 

As described in connection with screening non- 
perturbative effects become important when qo^T which 
corresponds to as^O.l. We shall treat the two limits sep- 
arately starting with the weakly interaction plasmas for 
which the gluon self energies, IlL,Ti^J■)^ are given by Eqs. 
(||) and (|9|). It is straight forward to calculate Q nu- 
merically and the result will be given below, but let us 
first make a simple analytical estimate. The main contri- 
bution to this integral can be obtained by including the 
leading terms in the self energies dsUy) 



QD: 



TlL{q,io) 



Uriq.uj) ~i-nql,. 



(B2) 
(B3) 



Thus we find for (Bl) 



Q 



= -n-a 

3^ 



qh 



ql 



■ In 1 



q^ 



TT 






7^ 9d 



2 
max 



(B4) 



Expanding in the limit qmax S> qo or equivalently for 
small as we obtain the leading orders up to a^ in the 
coupling constant 



Q 



In 



where 



1 / ,4 
qmin = qo exp <; - I 1 - In - 



> to, 



l.lSto 



The two terms in (B5) corresponding to Ing, 
Inq, 



(B5) 



(B6) 



and 



1.26gz 



(E 



because the additional terms in Hl^t lead to some adi 
tional screening besides the Debye screening and Land 
damping of (B2) and (B3). This effective cutoff is c 
termined by the screening only and is therefore the sai 
for gluon-gluon, quark-gluon and quark-quark scatterii 
Whereas qmin may serve as an effective "cutoff' of sni 
momentum transfers interactions, it is not a paramei 
put in by hand as discussed in pn]. Contrarily, it 
caused and determined by Debye and dynamical scret 
ing. 

If the transverse interactions are assumed to be Det 
screened like the longitudinal ones, i.e., 11^ = q^, th 
the result would have been qmin = qo exp(0.5) = 1.65(j 
This is because dynamical screening of Eq. (p3) is less 
fective than the Debye screening of (B2) and thus resu 
in a smaller qmin- 

It is convenient to express the results in terms of c 
In weakly interacting plasmas we find 



Q 



In 



2 



qn 



1.26qDj 

{qn 



-U 



■:/Tr 



Anasil + Nf/e) 



:0.1. (E 



The upper effective cutoff qmax is provided by t 
quark and gluon distribution functions as discussed 
connection with Eq. (|2^) and it therefore varies son 
what with particle type. Because Bose distribution fui 
tions emphasize smaller momenta than Fermi ones, g„ 
is larger for quarks. We find g^£ = 3.0T, qi^'J, = 3.8 
and qiriax = 4.8T for gluon-gluon, quark-gluon, a 
quark-quark scattering respectively. The lower effect: 
cutoff qmin is, however, the same for the three cases I 
cause it only depends on the screening in the gluon pre 
agator. Furthermore, we find that the extra terms in t 
matrix elements of (0), (0) and (0) besides the t^^ p; 
do not contribute much since they have varying signs a 
turn out to be partially cancelling. Thus the constai 
within the logarithms of ( p7|) and ( [40| ) just reflects t 
different qmax for gluon-gluon, gluon-quark and qua] 
quark scattering. 

Lacking screening of transverse interactions in t 
static limit, it has often been assumed that some mec 
anism like Debye screening might lead to screening 
transverse interactions as well, i.e. mpi = q^. Recent 
lattice gauge calculations of QCD plasmas have fou 
effective screening masses of order nipi ~ LIT near t 
phase transition point, Tc — ISOMeV. In both cases it 
thus assumed that 



Hl = nT/(l - Ai') = 



(E 



imin are often referred to as 
tributions in the literature |0 

A numerical evaluation of (pi 
(H) and (y) instead of ( B2 ) and ( B3|) gives a slightly larger 
value for the effective minimum momentum transfer 



hard" and "soft" con- in (Bl) which leads to 



with IiL,T given by Eqs. 



Q 



In 1 



q^ 



q^ 

^max 



m. 



pi 



mil + g^, 



(Bl 
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With qtfl 



o nrp (gq) 

(ss) 



o orp Jqq) 

o.o± , qrnax 



I.IT we find Q{q)^Vx/q. 

Aqq) 



4.8T and m 

(99) 



pi 



= 0.626, Qiqi^lL/qrmn) = 

0.819 and Q{qmax/qmin) = 1-024. These values enter the 
tts^LO.l expressions of Eqs. ( |3^ ) and (pl). 
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FIG. 1. Feynman diagram for gluon-gluon scattering in the 
t-channel. The hnes in the loops can be either quark or gluon 
propagators. 



FIG. 2. The collision geometry. For small momentum 
transfer, q -^ pi,p2, energy and momentum conservation re- 
quires cos 6i — cos 02 = i^/q- 



FIG. 3. The gluon viscosity for Nf — 3 assuming Debye 
and dynamical screening (dashed curve) , a constant screening 
mass TUpi = I.IT (dashed-dotted curve) and the minimal one 
(full curve). 



FIG. 4. The quark, gluon and total viscosities for Nf = 3. 
FIG. 5. The quark, gluon and total viscosities for Nf = 2. 



FIG. 6. The function f/x as determined by (A15). Also 



shown are the limits of (A16), (A17) and the simple ansatz 
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